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Abstract. We provide a derivation of the Foppl-von Karman equations for the shape 
of and stresses in an elastic plate with residual strains. These might arise from a range 
of causes: inhomogeneous growth, plastic deformation, swelling or shrinkage driven by 
solvent absorption. Our analysis gives rigorous bounds on the convergence of the three 
dimensional equations of elasticity to the low-dimensional description embodied in the 
plate-like description of laminae and thus justifies a recent formulation of the problem to 
the shape of growing leaves. It also formalizes a procedure that can be used to derive other 
low-dimensional descriptions of active materials. 

I. Introduction 

Laminae or leaf-like structures are thin, i.e. they have one dimension much smaller 
than the other two. They arise in science and technology in a variety of situations, from 
atomically thin graphene (thickness h ~ o(l)nm) with a lateral span of a few cm, to the 
earth's crust (h ~ 10km) which spans thousands of km laterally. On the everyday scale, 
recently there has been much activity on trying to understand the mechanics of these 
laminae when they are active or actuated, as in a growing leaf, a swelling or shrinking 
sheet of gel, a plastically strained sheet etc. In all these situations, the shape of the lamina 
arises as a consequence of the fact that inelastic effects associated with growth, swelling 
or shrinkage, plasticity, etc. results in a local and heterogeneous incompatibility of strains 
that leads to local elastic stresses. When combined with force balance, this naturally leads 
to the non-trivial shapes that are seen even in the absence of any external forces. A simple 
experiment suffices to make this point - when growing leaves or plastically strained ribbons 
[22J are cut in different directions to partially relieve the incompatible strains due to growth 
or plasticity, they relax to different shapes. 

Recently, these observations have lead to a quest for a theory that describes the coupling 
between residual strain that might arise from a multitude of causes to the ultimate shape of 
the object. Given that the deformations and strains involved are not necessarily small, this 
raises an age-old question of how to decompose the deformations into the elastic and inelas- 
tic parts. A possible approach is to consider the equilibrated shape of grown bodies and use 
a multiplicative decomposition of the deformation gradient that borrows from the theory 
of crystal plasticity which requires the notion of a reference configuration with respect to 
which all displacements are measured. This equilibrium approach leads to conceptual diffi- 
culties when matter is not conserved, as in growth processes, and suggests an incremental 
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approach embodied in an evolutionary rather than an equilibrium process, which also has 
antecedents in the theory of plasticity. In either case, we then require the knowledge of 
a constitutive law in addition to a characterization of the geometry of the body. For soft 
materials such as plant and animal tissue, a reasonable assumption is that insofar as the 
elastic response is concerned, the material is hyperelastic, while the inelastic deformations 
follow different laws depending on their origin. For example, swelling or shrinking gels may 
be described by a poroelastic theory that couples fluid flow to osmotic stress and deforma- 
tion, while growth in biological tissues arises from cell proliferation or cell shape change 
that have their own description which couple molecular and macroscopic processes. For 
example, at the molecular level, mutants responsible for differential cell proliferation [20J 
lead to a range of leaf shapes. At the macroscopic level, stresses induced by external loads 
lead to phenotypic plasticity in algal blades that switch between long, narrow blade-like 
shapes in rapid flow to broader undulating shapes in slow flow [T3] . 

Recent work has focused on some of these questions by using variants of thin plate 
theory to highlight the self-similar structures that form near the edge due to variations in a 
prescribed intrinsic metric of a surface that is asymptotically flat at infinity [T], and also on 
the case of a circular disk with edge-localized growth [3 [3], the shape of a long leaf [TS] etc. 
However, the theories used are not all identical and some of them arbitrarily ignore certain 
terms and boundary conditions without prior justification. This suggests that it might 
be useful to rigorously derive an asymptotic theory for the shape of a residually strained 
thin lamina to clarify the role of the assumptions used while shedding light on the errors 
associated with the use of the approximate theory that results. Recently, such rigorous 
derivations were carried out [7J [HI [151 E] i n the context of standard nonlinear elasticity 
for thin plates and shells. Further, in [16], a residually strained version of the Kirchhoff 
theory for plates [11] was rigorously derived under the assumption that the target metric 
is independant of thickness. 

In this paper, we carry out such a derivation under a different assumption on the asymp- 
totic behavior of the prescribed metric and show that the resulting equations are identical 
to those postulated to account for the effects of growth in elastic plates [18] and used to 
describe the shape of a long leaf. We limit ourselves to the case when a decomposition of 
the deformation gradient into an elastic and inelastic part can be carried out - this requires 
that it is possible to separate out a reference configuration, and is thus most relevant for 
the description of plant morphogenesis. Although our results are valid for thin laminae 
that might be residually strained by a variety of means, here we will limit ourselves to 
studying the role of incompatible growth strain in determining the shape of the laminae. 
In particular, we will only consider the one-way coupling of growth to shape and ignore the 
feedback from shape back to growth. 

We now proceed to give an overview of the main results in the paper. 

1. Geometry of incompatible strain. For a given mid-plate Q which is an open 
bounded and simply connected subset of R 2 , consider a sequence of 3d plates: 



Q h = Q x (-h/2, h/2), 0<h«l, 
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viewed as the reference configurations of thin elastic (and homogeneous) tissues. A typical 
point in Q h has the form (x',hx^) where x' G Vt and \x%\ < 1, and we shall make no 
distinction between points x' E fl and (x', 0) G Q h . 

Each Q h is assumed to undergo a growth process, whose instantaneous growth is described 
by a smooth tensor a h = [a^] : Q h — > M 3x3 , with the property: 

VxG£> deta h (x)>0. 

In general, a h may follow its own dynamical evolution. However, here we will focus only 
on the effect of the tensor a h on the effective elastic theory for the grown body. Here we 
will use a multiplicative decomposition of the deformation gradient that is similar to the 
one used in plasticity and also used in various growth formalisms, e.g. [21J, u : Q h — > UX 3 : 

V« = Fa h . 

The tensor F = Vu(a h )~ l corresponds to the elastic part of the deformation u, and accounts 
for the reorganization of the body Q h in response to the growth tensor a h . The above 
assumes that it is possible to differentiate a reference configuration with respect to which 
one might measure all relative displacements. As mentioned before, this is true for some 
but certainly not all growth processes. In particular, this is reasonable for botanical growth 
processes, but is unlikely for animal growth and remodeling processes that include large 
scale tissue flows and movements. 

2. Elastic energy associated with residual strains. The elastic energy of u is now 
a function of F only, and it is given by: 

(1.1) I^( u ) = \f W{F)dx = \( W{Vu{a h )- x ) dx, \/u G W 1,2 ({l h , M 3 ). 

h Jnh h J Q h 

The elastic energy density W : IR 3x3 — > M + is assumed to be compatible with the condi- 
tions of normalization and frame indifference (with respect to the special orthogonal group 
5*0(3) of proper rotations in IR 3 ): 

(*i) VF G M 3x3 Vi? G SO(3) W(R) = 0, W(RF) = W(F). 

Further, we shall require the nondegeneracy of W in the sense that: 

(*ii) 3c> VF G M 3x3 W(F) > c dist 2 (F, SO(3)), 

and also assume that W is C 2 regular in a neighborhood of 5*0(3). 

We note that by the polar decomposition theorem every a h can be uniquely written as 
a product R h a h , where R h G 50(3) and a h is symmetric positive definite. Hence, if W is 
moreover isotropic: 

(1.2) VF G R 3x3 VF G 50(3) W(FR) = W(F), 

then without loss of generality we can assume a h to be symmetric positive definite. 
As we shall see, it is instructive to study the following energy functional: 

(1.3) I%(u) = \( dist 2 (Vu(x)(a h (x))- 1 ,SO(3)) dx, 

h Jnh 

as the energy in (11. ip obeys a bound from below: Iyy > cIq. 
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3. Relation with the non-Euclidean elasticity. We shall here compare the above 
approach with the target metric formalism proposed in [5j and further developed in [16J. 
On each Q h we assume that we are given a smooth Riemannian metric g h = [g^]- The 

matrix fields g : Q h — > M 3x3 are therefore symmetric and strictly positive definite up to 
the boundary dQ h . Let \/~g^ be the unique symmetric positive definite square root of g h 
and define, for all x G Q h , the set: 

(1.4) J*(x) = ^Ry/g^{x); R G 50(3)} C M 3x3 . 

By the polar decomposition theorem, the necessary and sufficient condition for a deforma- 
tion u on Q h to be an orientation preserving realization of g : 

(Vu) T Vu = g h and detVw > a.e. in Q h 

is the following: 

Vu(x) G T h {x) a.e. in VL h . 
Motivated by this observation, we can replace the energy functional (jl.3p above by: 

(1.5) I%{u) = \( dist 2 (Vu(x),F h (x)) dx Vuelf 1 ' 2 ^,! 3 ), 

h Jnh 

measuring (in L 2 ) the pointwise deviation of u from being an orientation preserving re- 
alization of the given metric g . We note that Iq is comparable in magnitude with Jq. 
Indeed, the intrinsic metric of the material is transformed by the growth tensor a h to the 
target metric g h = (a h ) T a h and it is only the symmetric positive definite part of a h given by 
y/g^ = a/ (a h ) T a h which plays the role in determining the deformed shape of the material. 

Our main results can be divided in four major subcategories, presented below. 

4. Main results: Scaling analysis of thin non-Euclidean plates. Given a sequence 
of growth tensors a h , each close to Id and defined on fl h , the main objective is to analyse 
the behavior of the minimizers of the corresponding energies 1^ as h — > 0. Let us recall 
that, as proved in [IB] , the infimum: 

m fc = inf{i*(u); u G W x > 2 (Q h , M 3 )} 

must be strictly positive whenever the Riemann curvature tensor of the metric g h = (a h ) T a h 
does not vanish identically on Q h . This condition for g h , under suitable scaling properties, 
can be translated into a first order curvature condition (II. 6ft below. In a first step (Theorem 
II. ip we establish a lower bound on rrih in terms of a power law: 

rrih > ch 13 , 

for all values of (3 greater than a critical (3q in (11. 7p . This critical exponent depends on the 
asymptotic behavior of the perturbation a h — Id in terms of the thickness h. 

Under existence conditions for certain classes of isometries, it can be established that 
actually ~ h 130 . In other words, our analysis includes identification of the magnitude of 
the elastic energy of minimizers of 1^, in terms of the thicknesss h. 
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The following quantity measures the essential variation of the tensors a h : 
Var(a h ) = ||V ten (a h | n)IU«»(n) + \\d^a h \\ L ^ {Qh) 
together with their scaling in h: 



ui = sup < u: lim - — Var(a h ) — 1 . 



The symbol Vtan denotes taking derivatives di and d 2 in the in-plate directions e\ = 
(1,0, 0) T and e 2 = (0,1, 0) T . The derivative <9 3 is taken in the out-of-plate direction 
e 3 = (0, 0, 1) T . We will work under the following hypothesis: 

(♦i) H a/l ||L°°(n' 1 ) + ||( fl/l ) 1 \\L°°(n h ) < C, 

(4kii) wi > 

Theorem 1.1. Assume (4)- Assume that for some uo > 0, i/iere exzsfo £/ie feji: 

1 /-V2 

e g (x') = lim — f a fe (x', t) - Id dt in L 2 (tt, M 3x3 ). 
which moreover satisfies: 

(1.6) curFcurl (e s ) 2x2 ^ 
and t/iai wo < min{2u;i, a>i + 1}. Then, for every (3 with: 

(1.7) > A, = max{w + 2, 2w }, 

we /iawe: lim sup —7 inf 7q = +00. 

Above, we used the following notational convention which will be employed throughout 
the paper. For a matrix F, its n x m principle minor is denoted by F nxm . The superscript 
T refers to the transpose of a matrix or an operator. The operator curl T curl acts on 2 x 2 
square matrix fields F by taking first curl of each row (returning 2 scalars) and then taking 
curl of the resulting 2d vector, so that: curFcurlF = df 1 F 2 2 — dj 2 (F 12 + F 2 \) + d 22 F n . The 
symmetric part of a square matrix F is denoted by sym F = 1/2(F + F T ). In particular, 
we readily see that: curl T curl F = curl T curl(sym F). Physically this condition corresponds 
to the fact that the growth strain e g is incompatible, i.e. it is not uniquely integrable and 
thus is not derivable from an elastic deformation gradient. 

5. Main results: Compactness. By a compactness result we mean identification of 
the limit behavior of the minimizing sequence to 1^. More generally, this analysis can also 
be done for any sequence of deformations u h G W l,2 (Q h , M 3 ) whose energy J^ / (m /i ) scale like 
h^°. In the scaling regimes considered in this paper this compactness result has the following 
form: first, modulo rigid motions the deformations u h converge, up to a subsequence and 
in a suitable space, to the identity map on fl. Second, the suitably re-scaled displacement 
fields converge to elements of certain classes of Sobolev infinitesimal isometries. 

Note that no assumptions will be made on the special form of the deformations u h . From 
this point of view our analysis is Ansatz-free and the limiting behavior of minimizers is 
rigorously shown to depend only on the choice of the sequence a h . 
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We present the compactness result (Theorem 11.21 below) assuming the special form of 
the growth tensor (11.81) which corresonds to the von Karman model where /3q = 4. In 
this context, the out-of-plane displacement v will be rescaled by the thickness h, and the 
in-plane displacement w by h 2 . Again, these scalings are naturally imposed by the original 
choice of the growth tensor a . 

Theorem 1.2. Given two smooth matrix fields e g ,K g ■ — > M 3x3 , define the growth 
tensors as: 

(1.8) a h {x', x 3 ) = Id + h 2 e g (x') + hx 3 K g (x'). 

Assume that the energies of a sequence of deformations u h G W 1,2 (Q h , R 3 ) satisfy: 

(1.9) I^{u h ) < Ch\ 

where W fulfills (X) ■ Then there exist proper rotations R h G 5*0(3) and translations c h G R 3 
such that for the normalized deformations: 

y h {x', x 3 ) = (R h fu h (x', hx z ) - c h : ft 1 — > R 3 

the following holds. 

(i) y h {x\x 3 ) converge in W 1,2 ^ 1 , R 3 ) to x' . 

(ii) The scaled displacements: 

1 f 1/2 

(1.10) V h {x') = - f y h (x', t) - x' dt 

h J -i/2 

converge (up to a subsequence) in W 1,2 (fl, M. 3 ) to the vector field of the form (0, 0, v ) T , 
with the only non-zero out-of-plane scalar component: v G W 2 ' 2 (Q, M.). 
(hi) The scaled in-plane displacements h~ x V^ an converge (up to a subsequence) weakly in 
W 1,2 (Q,R 2 ) to an in-plane displacement field w G W 1,2 (Q, R 2 ). 

6. Main results: T-convergence. Heuristically, a sequence of functionals F n is said 
to T converge to a limit functional F if the the minimizers of F n , if converging, have a 
minimizer of F as a limit. More precisely, any T-convergence result involves a careful 
comparison of the values of the energies F n on seqeunces of deformations and the value 
of F on the limit deformation. Hence, it combines a lower and an upper bound estimate 
(which are called the T-liminf and the T-limsup inequalities). 

For the von Karman growth tensor studied in this article, these estimates are established 
for the sequence l/h A I^ v {u h ) and the limit energy value Z g (w,v) given in (II. lip below. The 
liminf inequality (Theorem ll.3p involves a lower bound on the energy of any sequence of 
deformations u h . The limsup part (Theorem IX .4j) establishes that for any pair of displace- 
ments (w,v) in suitable limit spaces, one can construct a sequence (14. ip of 3d deformations 
of thin plates Q h which approximately yield the energy I g (w,v). The form of such re- 
covery sequence delivers an insight on how to reconstruct the 3d deformations out of the 
data on the mid-plate Q. In particular, comparing the present von Karman growth model 
(14. ip with the classical model ([7], Section 6.1) we observe the novel warping effect in the 
non-tangential growth (see f l4.2p ). 
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Theorem 1.3. Assume ( II. gj) and Let the bound ( fl.^j) 6e satisfied by a sequence 

u h e ry 1 ' 2 (f2 h ) R 3 ) so t/iat £/ie convergences (i), (ii) and (Hi) of Theorem hold true. 
Then there holds: 



]immi^(u h )>l g (w,v), 



where: 



'1.11 



^g)2x2 I j 



Z fl (io, v) = ~ J Q 2 ^symVw + ^Vn <S> Vn - (sym e 9 ) 2x2 ^ 

+ ki L Q2 ( v2y + (sym s 

and the quadratic nondegenerate form Q 2 , acting on matrices F G R 2x2 is: 

(1.12) Q 2 (F) = min{Q 3 (^); F G R 3x3 , F 2x2 = F} and Q 3 {F) = D 2 W{ld)(F ® F). 

Theorem 1.4. Assume U.8\) and Then, for every w G W 1,2 (Q, R 3 ) and every 

n G W^ 2 ' 2 (f2, R) ; there exists a sequence of deformations u h G W X,2 {S \ R 3 ) snca iaai iae 
following holds: 

(i) T/ie sequence y h (x',x 3 ) = u h (x', hx^) converge in VF 1 ' 2 (0 1 , R 3 ) to x'. 

(ii) F'V) = /T 1 -/- (« h (a/, t) - x') dt conwerae in W 1 ' 2 ^, R 3 ) to (0, 0, v) T . 

J-h/2 

(iii) h^V^ converge in W /1,2 (f2,R 2 ) to w. 

(iv) Recalling the definition ( CQZj ) one /jas: 

£S^Sr(« fc )=V«'.«)- 

The main consequence of the T-convergence result is the following: If u h is a minimizing 
sequence for Iyy, and if if and v are the respective limiting in-plane and out-of-plane dis- 
placements corresponding to u h , then (w, u) will be a minimizer of the von Karman growth 
functional T g (Corollary 11.51 below) . 

Another direct corollary is the identification of the lower bound of for the von Karman 
growth (11.8j) under an appropriate curvature condition. Note that the assumptions of 
Theorem 11.11 do not hold, since in the present case u = 2u>i = U\ + 1 = 2. However, if we 
replace (11 .6p by either of the following two conditions: 

(1.13) curl((sym k 9 ) 2x2 ) ^ 0, 



or: 



(1.14) curFcurl (e 9 ) 2x2 + det((sym /%) 2x2 ) ^ 0, 

we indeed obtain that inf 1^ > ch A with c > 0. The above conditions guarantee that the 
highest order terms in the expansion of the Riemann curvature tensor components #1213, 
-R2321 and -R1212 of 9 — {a h ) T a h do not vanish. Also, either of them implies that inf T g > 
(see Lemma [6. ip . which combined with Theorem 11.31 yields the lower bound on inf l\. 

Corollary 1.5. Assume U.8\) and (Jfr). Then: 
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(i) There exist uniform constants C,c>0 such that for every h there holds: 

(1.15) c<^MI^<C. 

If moreover $1.13\) or ( 1.14 ) hold then one may have c > 0. 

(ii) There exists at least one minimizing sequence u h G W /1 ' 2 (fi /l ,M 3 for 1^: 

(1.16) £ j(ii*(„»)-^mfi*)=0. 

For any such sequence the convergences (i), (ii) and (Hi) of Theorem \1.2\ hold and 
the limit (w,v) is a minimizer ofT g . 

(iii) For any minimizer (w,v) of T g , there exists a minimizing sequence u h , satisfying 
M.1Q) together with (i), (ii), (iii) and (iv) of Theorem \1.4 



7. Main results: Euler-Lagrange equations for the limit theory. When W is 
frame invariant as in (Jfri) and isotropic (II. 2p . one can see [7] that the quadratic forms of 
(j!.12p are given explicitely as: 

(1.17) Q 3 (F) = 2/i|sym F| 2 + A|tr F| 2 , Q 2 (F 2x2 ) = 2/i|sym F 2x2 | 2 + -^-|tr F 2x2 | 2 , 



2fi + X 

for all F G IR 3x3 , where tr stands for the trace of a quadratic matrix, and fi and A are the 
Lame constants, satisfying: /i>0, 3A + /i>0. 

We will show that under these conditions, the Euler-Lagrange equations (15.31) . (15.81) of T g 
in (11.111) are equivalent, under a change of variables which replaces the in-plane displacement 
w by the Airy stress potential $, to the system of von Karman-like equations introduced 
recently by Liang and Mahadevan in [18J: 

A 2 ^ = S{KQ + \g) 



(1.18) 



BA 2 v = [v, $] - BVt 



fl ! 



under the corresponding boundary conditions (I5.12p (I5.13P (I5.14p . where: 

S= Young's modulus _ M(3A + /i) ^ Kg = Gaussian curvature =-[t>,t>], 

A + jj, 2 

S A 
B = bending stiffness = — v = Poisson's ratio 



12(1 -u 2 Y 2(A + /i)' 

X g = curFcurl (e 9 ) 2x2 = 9 22 (e s )n + d u (e g ) 2 2 - <9i 2 ((e 9 )i 2 + (e s ) 2 ij, 

Q g = div T div (^K g ) 2x2 + v cof (k 9 ) 2x2 ^ 

= (9n f(/« g )ll + v(Kg)22 \ + 8 2 2 ((Kg) 2 2 + v(Kg)ll) + (1 ~ v)d X2 {[^g)\2 + (%) 



21 



In brief, we may summarize our results as follows. Under the special form of the growth 
tensor given by (II .8p . we have: 

(i) rrih = inf 1^ scales like h 4 , 
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(ii) there exists a sequence of deformations u h such that Iy^(u h ) — rrih = o{h A ) 1 and 
for which the rescaled in-plane and out-of-plane displacements converge to a limit 
(w,v), 

(hi) (w, v) minimizes the von Karman growth functional X g in ( II. lip and hence satisfies 
the system ( I1.18P with the corresponding free boundary conditions ( I5.12p ( 15 . 1 3[) 
( I5.14p . when expressed in terms of the Airy stress potential. 

8. Approximating low energy deformations. A crucial step in obtaining the above 
results is an approximation theorem (Theorem 11.61 below). The underlying idea is that we 
can control the oscillations of the deformation gradient Vu h in boxes of diameter propor- 
tional to the thickness h. The geometric rigidty estimate of Friesecke, James and Miiller 
[6] and its generalization to our setting is the basic tool in this step. 

Theorem 1.6. Assume (4)- Let u h E W h2 (Q h , M 3 ) satisfy: 

hm ^J V) = 0. 
h^o h 2 

Then there exist matrix fields R h E W 1>2 (£1, R 3x3 ) ; such that R h (x') E SO (3) for a.a. x' E tt 
and: 

1 ; \Vu h (x) - R h {x')a h {x)\ 2 dx<C (l^(u h ) + h 2 Var 2 (a h )) , 



h 



n h 



\VR h \ 2 < Oh' 2 (l£(u h ) + h 2 Var 2 {a h )) , 
where the constant C is independent of h. 

2. Scaling analysis: a proof of Theorem 11.11 

1. Take f3 > max{w + 2,2w } and assume, by contradiction, that for some sequence 
u h E W^'^^M 3 ) there holds: 

(2.1) ^loV) < C. 

Since (3 > 2, in virtue of Theorem 11.61 there exists the rotation-valued matrix fields R h E 
W 1,2 (fl, SO(3)) approximating appropriately V'U /l (a' l )~ 1 . Observe that: 

(2.2) dist 2 (-f R h ,S 0(3)) < 1 \R h {x)-R h {x Q )\ 2 dx<C [\VR h \ 2 ^0 zsh^O, 

Jn Jn Jn 

by the second estimate in Theorem 11.61 and by (4&). We may hence, for small h, define the 
averaged rotations R h E 50(3) by R h = F so(3) f Q R h . 

Define now two fields: V h E W 1,2 (fi,]R 3 ) and A h E W^ 2 {Q, M 3x3 ): 

1 rh/2 

V\x') = -—f (R h ) T u h (x', t) - x' dt, 
h"° J- h ,2 

A h (x') = \(R h ) T R h (x') ^ a h (x',t) dt - Id ) . 



-h/2 
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Observe that: 



n r rh/2 

W h - A% x2 \\ 2 L2{a) < — ^ j ^R h (x')a h 3x2 (x',t)-V tan u h (x',t) dt 



dx' 



< 



C 



\Vu h (x) - R h (x')a h (x)\ 2 dx 



< C ( T^oV) + —^Var 2 (a h )] ^ as h 0, 

by Theorem 11.61 and since 2u < (3 and 2oj — 2 < 2u\ . 
2. Notice that: 

A h (x') = (R h ) T R h (x')e g (x') + ((R h ) T R h (x') - Id) 



- (R h ) T R h (x') e g (x') 



h/2 



h/2 



a h (x',t) - Id dt 



Clearly, the third term above converges in L 2 (Q) to 0, by the definition (11. ip . The first 
term converges to e g , as by (12. 2p and Theorem 11.61 



(2.3) 



\(R h ) T R h -W <C I \R h -R 



bh\2 



< C ( [ \R h -i R h \ 2 + dist 2 (-/ R h , SO (3))) <C f \VR h \ 2 ^ as h -> 0. 
\Jn Jn Jq J Jn 

To deal with the second term, in the expansion of A h , recall that R h ,R h (x') G 50(3) 
and so: 

1 -sym ((R h fR h - Id) = {{R h ) T R h - \df ((R h ) T R h - Id) . 



Therefore, reasoning as in in (12. 3p : 
1 



^\\sjm((RrR h -ld)\\ L2{Q) < ±\\(R h ) T R h -ld\ 

C 



2 

L 4 (f2) 



c 



<< ~ II J?' 1 r?ft- 1 1 2 ^ 



< 



C (-±—I%(u h ) + -^-Var 



VjR h|2 

a ft ) ] -»■ as /i 0, 



n 



by Theorem 11.61 and since Co> + 2 < j3 and cj < 2oj\. 
3. Summarizing, the previous step yields: 



sym e g = lim sym A h in Z/(fi, M Jxd ), 



which by O) implies: 



(sym e 9 ) 2X 2 = lim sym V^' 1 in L 2 (£l, \ 

h— s-0 



>2x2\ 



Consequently, (sym e 9 ) 2x2 = symVV, for some V G W /1,2 (fi, M 3 ), and hence there must be: 
curl T curl (e 9 ) 2x2 = curl T curl (sym e 9 ) 2x2 = in fi. This brings a contradiction with (I1.6P 
and hence ends the proof. ■ 
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3. Lower bound for the Von Karman scaling: a proof of Theorem 11.21 and 

Theorem 11.31 

Consider a specific form of the growth tensor: 

a h ( x ', x 3 ) = Id + h~<e g (x') + h 8 x 3 K g (x') 

with exponents 7, 9 > and the smooth perturbation moments: e g , K g : Q — > M 3x3 . One 
easily sees that (4) holds here with: 

Var(a h ) = h 1 1| Ve 9 1| l°° + h 9 \\ n g \\ LOO , coi = min{7, 9}, co = 7. 

Also, the result in Theorem 11.11 holds provided that curl T curl (sym e g ) ta n ^ and: 

(3.1) 7 < min{0 + 1,29} and /3 > max{ 7 + 2, 2 7 }. 

In what follows, we shall work with exponents j = 2,9 = I, ft = 4 which are critical for 
both inequalities in ( 13. ip . 

Proof of Theorem 11.21 

1. Let R h G W 1,2 (fi, 50(3)) be the matrix fields as in Theorem [L6J 

(3.2) \ [ \Vu h - R h a h \ 2 < Ch 4 , [ \VR h \ 2 < Ch 2 . 



h 



Define the averaged rotations: R h = F S o(3) f n R h - These projections of j R h onto 50(3) 
are well defined for small h in virtue of: 

dist 2 (-/ R h , 5-0(3)) < / \R h (x) - R h (x )\ 2 < C [ \VR h \ 2 < Ch 2 . 
Jn Jq Jn 



Further: 



(3.3) J \R h - R h \ 2 <c[j \R h - J R h \ 2 + dist 2 (J R h ,S0(3))^j < Ch 2 
Let now: 

(3.4) R h = P 50( 3) / {R h ) T Vu h . 

The above projection is well defined for small h, because dist 2 (f nh (R h ) T Vu h , 50(3)) is 
bounded by: 

I -f (R h ) T Vu h - Id| 2 < C i \Vu h -R h \ 2 
Jn h Jn h 

<c(i \Vu h - R h a h \ 2 + £ \a h -\d\ 2 +l \R h - R h \ 2 ) < Ch 2 , 
\Jn h Jn h Jn h ' 

where we used f)3.2p and (13. 3p . Consequently, we also obtain: 

(3.6) \R h - Id| 2 < 0|skew i {R h ) T Vu h \ 2 <C\1 (R h ) T Wu h - Id| 2 < Ch 2 . 

The first inequality above follows by noticing that for any matrix F close to Id there holds: 
P 50 (3)(sym F) = Id, and hence: |P SO ( 3 )F - Id| < C\F — sym F < 0|skew F\. 



(3.5) 
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2. We may now define: 

(3.7) R h = R h R h . 
By d33D, AMD and (Q it follows that: 

(3.8) f \R h - R h \ 2 <Ch 2 and lim(£ A ) T fl & = Id in W 1 ' 2 ^, M 3x3 ). 

Let c ft G K 3 be vectors such that for the rescaled averaged displacement V h defined as in 

,/2 
/i/2 l 



(ILTO]) : V*(a/) = /i- 1 ^i 2 JR h ) T u h (x',t) -c h -x' dt, there holds 



(3.9) / V h = 0, skew / W ft = 0. 

in Jn 

The second statement in ( 13.91) follows by noticing that, for a matrix F sufficiently close 
5*0(3), its projection R = ¥so(3)F is coincides with the unique rotation appearing in the 
polar decomposition of F, that is: F = RU with skew U = 0. Therefore and in view of 
(E3D and (EZD we obtain that (R h ) T j Qh Vu h = (R h ) T f Qh (R h ) T Vu h is symmetric. Hence: 
skew f n W h = /i _1 skew j- Qh (R h ) T 'Vu h = 0. In particular, we see as well that (13.71) coincides 
with: 

R h = ¥so(z)i V«" 
To obtain (i) we use (13. 7ft . ( 13.51) and ( 13.61) : 



n* 



(3.10) 



||(V/-Id) 3x2 ||i w < i^J(^) T V^-Idp 



,,l (a;) — x' 



<c(^j \(R h ) T Vu h - Id| 2 dx + \R h - ld\ 2 ^j < Ch 2 
and notice that by (13. 2|) : 

||9 3 /||i w < ch [ \vu h \ 2 <ch 2 , 

Jn h 

which implies convergence of y h by means of the Poincare inequality f Q1 y h 
h j n V h = by ( 13. 9p . Notice also that (I3.10p implies the weak convergence (up to a 
subsequence) in H/ 1,2 (n,R 3 ) of V h . 

3. Consider the matrix fields A h e W 1,2 (ft, M 3x3 ): 

i rV2 

A V) = - f (R h ) T R h (x')a h (x', t) - Id dt 

" J-h/2 

= h(R h ) T R h (x')e g (x') + i ({R h ) T R h (x') - Id) . 
By (13. 8p and (13.21) . clearly: || A^Hvpwn) < C and so, up to a subsequence: 

limA' 1 = A and lim - ((R h ) T R h — ld)=A 
(3.11) fc->-o fc->o h KV J ' 

weakly in W 1,2 (fi, M 3x3 ) and (strongly) in L q (Q, M 3x3 ) Vg > 1. 
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Also, using (13 .8p and (13. 2 j) again: 

/^Hsyin^)^ - Id)|| i2(Q) = 2h- 1 \\{{R h ) T R h - \df{{R h ) T R h - Id)|| L2(c) 

< GhT x || {R h ) T R h - ldf L A m < Ch- 1 ^ 1 - R h \\ln,2(n) < Ch. 
Above, we used a straightforward observation that: 

(3.12) (R-ldf(R-ld) = -2sym(i2 - Id) VReSO(3). 
Therefore, the limiting matrix field A has skew values: 

(3.13) sym A = lim sym A h = 0. 

h— i>0 

Further, by (13 . 1 2[) we observe that: 

^sym A h = sym ((R h ) T R h e q ) - ~ ((R h ) T R h (x') - Id) T ((R h ) T R h (x') - Id) 
h 2h 2 

Hence, by (1EB|) . ( 13~TTD and (13TT3D : 

(3.14) lim ^symA' 1 = sym e 9 - -A T A = sym e 9 + -A 2 in L q (tt, M 3x3 ) Vg > 1. 
h->o h 2 2 

4. Regarding convergence of V h , we have: 

(3.15) = A h 3x2 (x') + -(R h ) T f R h (x')a h 3x2 (x',t) - V tan u h (x',t) dt. 

11 J-h/2 

Also: 

2 

dr' 

' -A/2 



c r r^ 2 

W* - A h 3x2 \\ 2 Lm <- j R h (x')a% x2 (x',t) - V tan u h (x',t) dt 

n Jfl J-h/2 



< £ / \Vu h (x) - R h (x')a h (xW d.r < Clr 



h 3 



and hence by (13.1 ip W h converges in L 2 (ft,M 3x2 ) to A. Consequently, by (l3~9|) : 
(3.16) lim V h = V in W 1,2 (Q, R 3 ), 7 G W 2 ' 2 (Q, M 3 ) and W = A 3x2 . 

h— >0 

Use now (13 . 1 3[) to conclude that symV(VJ an ) = and so by Korn's inequality V tan must be 
constant, hence in view of (I3.9p . This ends the proof of (ii). 

To deduce (iii), divide both sides of ( 13 . 1 5|) by h and pass to the limit with its symmetric 
part. By (13. 14j) . the first bound in (13 .2p and (13. 9p we conclude that: 

Wh-'VUw^n) < C||V(/r-V^)|b(n) = C\\sjmV(h- l V t h an )\\ LHn) < C, 
which proves the claim. ■ 

Proof of Theorem 11.31 

1. Define the rescaled strains G h G L 2 (fi 1 ,R 3x3 ) by: 

G h {x', x 3 ) = ^ ((R h (x')) T Vu h (x', hx 3 )a h (x', hx^' 1 - Id) . 
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Clearly, by (13 .2h \\G h \\yi(pi') < C and hence, up to a subsequence: 

(3.17) lim G h = G weakly in L 2 (n 1 ,R 3x3 ). 

h— >0 

We shall now derive a property of the limiting strain G. Observe first that: 

(3.18) lim -^(d 3 y h - he 3 ) = Ae 3 in L 2 {Q\ R 3 ), 
where e 3 = (0, 0, 1) T . This is because: 

±(d 3 y h (x)-he 3 ) = ^(R h ) T Vu h (x',hx 3 )-Idy 3 

^{R h ) T (Vu h {x', hx 3 ) - R h (x')a h (x', kc 3 ))e 3 

+ hR h ) T R h {x')(a h {x',hx 3 ) - Id)e 3 + U{R h ) T R h {x') - Id)e 3 , 

where the first term in the right hand side converges to in L 2 (fi x ) by ( 13. 2p . the second 
term to in L°°(fi 1 ), and the last term to Ae 3 in L 2 (f2) by (pTTT]) . 

For each small s > define now, with a small abuse of notation, the sequence of functions 
/^G^n 1 ,! 3 ): 

(3.19) f s > h (x) = ~(y h (x + se 3 ) - y h (x) - hse 3 ). 

Clearly f s ' h (x) = ± f° y h {x + te 3 )- he 3 dt, and so by fl3~T8|) : 

lim/ s,/l = Ae3 in L 2 {Q\R 3 ). 

ft— >o 

Also d 3 f s,h (x) = - s -^{d 3 y h (x + se 3 ) — d 3 y h (x)) so again by (I3.18P : 

\imd 3 f s ' h = in L 2 (n\R 3 ). 

Further, for any a = 1, 2 we have: 
1 1 



<9 a r fc (x) = -^-(R h ) T {Wu\x' } hx 3 + hs) - Vu h (x', hx 3 )y c 



1 (R h ) T R h (x') I G h (x', x 3 + s)a h (x', hx 3 + hs) - G h (x', x 3 )a h (x', hx 3 ) 



s 



+ -^{a h {x' ,hx 3 + hs) - a h (x\hx 3 )yje ai 

which, in view of H3. 1 71) and (13. 8p yields the weak convergence in L 2 (Q 1 ,R 3x2 ) of: 

(3.20) lim d a f s ' h (x) = - (g{x , x 3 + s) — G{x , x 3 ) j e a + K g (x')e a . 

Consequently, we see that f s,h converges weakly in W l,2 (Q, R 3 ) to Ae 3 . Hence, the left 
hand side in (I3.20p equals d a (Ae 3 ) and so: 

(3.21) G(x) 3x2 = Gq(x') 3x 2 + X 3 Gi(x') 3x 2, 



ELASTIC PLATES WITH RESIDUAL STRAIN 15 
for some G G L 2 (fi,R 3x3 ) where: 

(3.22) G 1 {x') = V(A(x')e 3 ) - k 9 {x'). 

2. Divide now both sides of ( 13.1 5p by h and pass to the weak limit in L 2 (Q,M. 3x2 ) 
with its symmetric part. Since lim h~ 1 symVV h = symVw by Theorem 11.21 (iii), and 
lim /i _1 sym A\ x2 = (sym e 9 ) 2x2 + l(A 2 ) 2x2 by (|3.14p, and: 



limsym [ (R h ) T R h (x') f ' G h (x' ,t)a h (x' , ht) dt] = G (x') 3>2 . 

^ 1 ' " ' 3 x2 



-1/2 
-1/2 

by dSUTD, (|3~2T|) and flU]), we obtain: 



1 



(3.23) symVw - (sym e s ) 2x2 ~ ^(^4 2 )2x2 = -G (x') 2x2 . 

3. We shall now prove the bound of the Theorem 11.31 First, Taylor expanding the 
function W(F) close to F = Id, and recalling (I1.12p we obtain: 

^(Vm^x)^^)- 1 ) = ^w(R h {x) T Vu h {x)a h (x)- A j 

= ^W(ld + h 2 G h (x)) = ^Q 3 (G h (x)) + h 2 0(\G h (x)\ 3 ). 

Consider now sets Qh = {x G fi 1 ; h\G h (x', h 3 )\ < 1}. Clearly xn h converges to 1 in L 1 (f2 1 ), 
with h — > 0, as hG h converges to pointwise a.e. by (13.21) . We get: 

liminf lrlk(u h ) > liminf i / X n h W(Vu h (x', hx 3 )a h (x', hx^)' 1 ) dx 

(3.24) = liminf Q Q 3 ( X n h G h ) + o(l) |Gf ) 

> - [ Qa(sym G(x)) dx, 

where the last inequality follows by ( 13. 2p guaranteeing convergence to Oof the termo(l) / \G h \ 2 , 
and by the fact that Xn h G h converges weakly to G in L 2 (f2 1 , M 3x3 ) (see (I3.17P ) in view of 
Q 3 being positive definite on and depending only on the symmetric part of its argument. 
Further: 

If If 

Q 3 (sym G)>- Q 2 (sym G 2x2 (x)) dx 



2 



= \ I ^ 2 ( sym ^0(^)2x2 + a; 3 sym G 1 (x') 2X 2^j dx 

(3.25) =-/ Q 2 (sym l7 (x / ) 2 x2) + r / x 2 Q 2 (sym G x {x') 2x2 ) 

1 Jn 1 1 Jn 1 

= ^ J Q2 ^sym Vw - (sym e 9 ) 2x2 - ^(^ 2 )2x2^j 

+ ^ J ^(sym (VAe 3 ) 2x2 - (sym /%) 2x2 ). 



16 



MARTA LEWICKA, L. MAHADEVAN AND REZA PAKZAD 



by 11531]) . (ESD and 1)522]) . Now, in view of Theorem O (ii) and ( 13361) one easily sees 
that: 



{A 2 ) 2x2 = -Vf ® Vv and 
which yields the claim by ( 13. 24ft and (j3.25p . 



VAe-- 



-W 



2x2 



4. Recovery sequence: a proof of Theorem 11.41 



For any F G 



d2x2 



by (F)* e 



d3x3 



we denote the matrix for which (F) 



2x2 



F and 



(F)* 3 = (F)* 3i = 0, % = 1..3. Recalling (fl~T2l . let c(F) e R 3 be the unique vector so that: 

Qz(F) = Q 3 ((i 71 )* + sym(c®e3)). 

The mapping c : ffi^m — ^ ^ 3 * s well-defined and linear, as Q 3 is a quadratic form, positive 
definite on the space of symmetric matrices. Also, for all F e R 3x3 , by 1(F) we denote the 
unique vector in R 3 , linearly depending on F, for which: 

sym(F - (F 2x2 )*) = sym(/(F) ® e 3 ). 

1 . Let the in-plane displacement u> and the out-of-plane displacement v be as in Theorem 
11.41 We first prove the result under the additional assumption of w and v being smooth up 
to the boundary. Define the recovery sequence: 



(4.1) u h (x',x 3 ) 



X' 





+ 



h 2 w(x') 
hv(x r ) 



+ x 3 



-hVv(x') 
1 



+ h 2 x 3 d°(x') + -hxl(f(u ) 



(4.2 



where the smooth warping fields d°, d 1 : Q — > R 3 are given by: 

d° = l(e g ) - ^|Vw| 2 e 3 + c^symVto - ^Vv <S> Vv - (sym 63)2x2), 
d l = 1(k 9 ) +c(~ V 2 v - (sym K g ) 2 x2\ • 

The convergence statements in (i), (ii) and (iii) of Theorem 1 1.41 are verified by a straightfor- 
ward calculation. In order to establish (iv) we need to estimate the energy of the sequence 
u . Calculating the deformation gradient we first obtain: 

Vu h = ld + h 2 {Vw)* + hA-hx 3 {V 2 v)* + h 2 [x 3 Vd° d°]+h[\x 2 Vd 1 x 3 d l ] , 

where the skew-symmetric matrx field A is given as: 



A 



-(Vv) T 
Vv 



We shall use an auxiliary 5*0(3)-valued matrix field R h = e hA . Clearly: R h = Id + 
hA + ^A 2 + 0(h 3 ) and (R h ) T = Id - hA + % A 2 + 0(h 3 ). Also, recall that: (a 71 )" 1 = 
Id — h 2 e g — hx 3 n g + 0(h 3 ). We hence obtain: 

{R h f{Vu h ){a h )- 1 = ld+h 2 ({Vw)*-l-A 2 -e g +d%e 3 )+hx 3 (-{V 2 vy- Kg +d 1 ®e 3 ) + 0{h 3 ). 
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Recalling now the definition of the quadratic form: 0.3(F) = D 2 W (Id) (F<g>F) = Q 3 (sym.F), 
Taylor expanding the energy density W around the identity, and taking into account the 
uniform boundedness of all the involved functions and their derivatives we get: 



e 3 



Ih r (u h ) = \( W(Vu h (a h )- 1 ) = l f w((R h ) T (Vu h )(a h ^ { 
h Jn h h Jn h v 

= y J Qs (sym((V™)* - \a 2 - e g + d° 

+ 24 / 2a (sym( - (V 2 w)* - K g + ® e 3 )) + 0(/i 5 ). 
Note that A 2 = (Vw <g> Vw)* - |Vw| 2 (e 3 <g> e 3 ). Therefore: 

sym ^(Vw)* - ^A 2 - e g + d° <g> e 3 ^j = ^symVw + ^Vw <g) - (sym e 9 ) 2X 2^ 

+ sym (^(d - l eg + i|Vw| 2 e 3 ) ® e 3 ^) , 

sym (-(V 2 w)* - Kg + d 1 <g> e 3 ) = (-V 2 w - (sym k 9 ) 2X 2)* + sym ((d 1 - l Kg ) <g> e 3 ) . 
In view of (14. 2[) it follows that: 

(4.3) ^(u h )=l g (w,v) + 0(h), 

which proves the desired limit (iv) for smooth displacements w,v. 

2. In order to carry out the analysis for w G W 1,2 (Q, M 2 ) and v G W 2,2 (Q,M.), it suffices 
to suitably approximate them in their respective norms by smooth sequences Wh and tv 
Define the sequence u h as in (14. ip using Wh and Vh instead of w and v for each h. The 
error O(h) in the final estimate (14. 3p . contains now an additional term hC(wh,Vh) where 
the quantity C{wh,Vh) depends only on the higher norms of Wh and Vh- This quantity can 
always be controlled by a uniform constant, by slowing down the rate of convergence of the 
sequences Wh and Vh- ■ 

5. EULER-LAGRANGE EQUATIONS OF THE FUNCTIONAL X g : THE DERIVATION OF 
(11.181) AND THE FREE BOUNDARY CONDITIONS 

Assume that (w,v) is a local minimizer of X g in (11. lip with Q 2 as in ( I1.17jl . 

1. Consider a variation <ft G C^°(Q, IR 2 ) in w. That is, for all small (positive or negative) 

e: 

l g (w + E(f), v) - X g (w, v) > 0. 
Collecting terms of order e, we obtain that: 

(5.1) f ( 2/i(sym V0) : (sym Vuj + ^) + 2/i \ (div0) (div w + tr ^) J = 0, 
where \1> G L 2 (fi,IR 2x2 ) is the vector field: 

(5.2) ^ = ivw (8) Vw - (sym e g ) 2x2 . 
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After integrating ( 15.1ft by parts and recalling the fundamental theorem of calculus of vari- 
ations, we obtain: 

(5.3) divM = 0, with M = 2/x(sym Vw + ^)+ 2 ^ (div w + tr #)Id. 

2/i -|- A 

Above, the divergence of a symmetric matrix field M is taken row-wise. Consequently, 
the z-th row of M (z = 1,2) can be written as V ± ip l for some scalar fields (ip 1 ,^ 2 )- The 
symmetric matrix field cof M has hence the form V(^ 2 , —ip l ) T , which implies that: 

cof M = V 2 $ for some $ G W 2 ' 2 (tt, R). 

Recall that for a matrix M G R nxn , cof M denotes the matrix of cofactors of M, that is 
(cof M)ij = detMij, where M tj G R^M™- 1 ) is obtained from M by deleting its 

ith row and jth column. 

From this discussion we see that (15. 3p is equivalent with: 

(5.4) M = cof V 2 $. 

In classical elasticity, the scalar field $ is called the Airy stress potential. 
2. We shall now need the following result: 

Lemma 5.1. Let a, (3 G R be such that: a ^ and a + 2(3 ^ 0. Then the following 
conditions are equivalent, for any matrix field F G L 2 (Q, R 2x2 ): 

(i) F = a symVw + (3(div u>)Id ; for some w G W 1,2 (Q, R 2 ), 

(ii) curl T curl B — — A(trF) = ; in the sense of distributions. 

a + 2/3 

Proof. We shall use the following easily obtained formulas: 

curl T curl (7M) = A7 = div T div (7W), 

valid for any scalar field 7 on O. 

To prove the implication (i)=^(ii), note that by (i): tr F = (a + 2(3) div w. Thus: 

curl T curl F = (3 curl T curl ^(div w)ldj = (3 A(div w) = — +2(3 / ^ r ^)' 
To prove the reverse implication (ii)=^(i) observe that by (ii): 

curl T curl F — curl T curl ( (tr -FUd^ = 0. 

a + 2(3 V ' J 

Thus: 

(5.5) F — -(tr F)Id = a sym Vw, 

a + 2(3 

for some vector field w. In particular: 

div w = tr (sym Vw) = — ( tr F tr F\ = -tr F. 

K 3 ' a\ a + 2(3 J a + 2(3 

Together with (15. 5ft the above implies (ii). ■ 
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We now use Lemma [5.11 with: 

F = cof V 2 $ - [2fiV + 2/iA - (tr mid] , a = 2u, (3 = 2uX/(2fj l + A). 
V 2 ix + A / 

Clearly, the condition (i) is equivalent to ( 15. 4p and hence (15. 4p is further equivalent to (ii), 
which after recalling (I5.2p takes the form: 

curFcurl (cof V 2 $) - - — ^— — A(tr cof V 2 $) = /x curFcurl (vv g) Vv - 2(sym e 9 ) 2x2 Y 
V / 2/x -\- oA V / 

Since curl T curl (sym e g ) 2x2 = curl T curl (e 9 ) 2x2 , an d curl T curl (Vv £g> Vv) = —2 det V 2 v 
and both curl T curl (cof V 2 $), and A(tr cof V 2 $) equal A 2 $, we obtain: 

2(/x + A) 



A 2 $ = -2fx(detV 2 v + curl T curl(e 9 ) 2x2 



2/i + 3A 
or equivalent ly: 

(5.6) A 2 $ = s(detV 2 v + curl T curl(e 9 ) 2x2 Y 



3. Consider now a variation tp G C^°(fl, R) in v, so that for all small positive and negative 

e: 

l g (w, v + Sip) - l g (w, v) > 0. 
Collecting terms of order e and calling: 

(5.7) 'I' = V 2 v + (sym k 3 ) 2x2 , 
we obtain that: 

(5.8) J ((Vcp <g> Vv) : M + BV 2 p : (^ + u cof = 0, 
where we used the following identity, valid for any F G M 2 ^: 

2//F + 2/xA (tr F)Id = 12£(F + i/ cof F). 

2yU -|- A 

By (15.41) the first term in the integrand of ( 15. 8ft equals V<^ • ((cof V 2 $)Vt>). Integrate by 
parts in ( 15. 8ft and use fundamental theorem of calculus of variations to obtain: 

(5.9) - div((cof V 2 $)Vt/) + 5 div T div + v cof = 0. 
Use now the following formulas: 

div((cof V 2 $)W) = (cof V 2 $) : V 2 v, div T div (V 2 v) = A 2 v, 

and remember that divcof of a gradient of a vector field vanishes, to find the following 
equivalent form of ( 15. 9ft : 

(5.10) (cof V 2 $) : V 2 v = BA 2 v + B div T div(jsym /%) 2x2 + v cof (sym /%) 2x2 ) 
Recalling the definition of Airy' bracket: 

= V 2 w : (cof V 2 $), 

so that det V 2 t> = l/2[v,v], we see that (15. 6ft and ( 15.1 Oft give the system (11. 18ft . 
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4. We will now derive the natural (free) boundary conditions satisfied by the minimizers 
(w,v) of X g in (II. lip . The analysis is equivalent as the above, for the variations <fi and <p 
which do not vanish on the boundary of dQ. 

Integrating (15. ip by parts and taking into account (15.31) . we obtain: 

(5.11) Mn = ondQ, 

where n denotes the normal to dQ. Hence by (15 .4p : (cof V 2 $)n = 0, which is equivalent 
to: <9 T V$ = 0, for the tangent vector field r to dfl. Therefore: 

V$ = const. on dQ. 

Since $ is detemined up to affine functions, we may assume that $(xo) and V$(a;o) vanish 
at a given point x £ 9Q. We obtain hence the first set of boundary conditions, for $: 

(5.12) $ = <9 rt $ = on dfl. 

To deduce the boundary conditions for the out-of-plane displacement v we use (15. 8p which 
is again valid for all ip G C°°(fi, E). Integrating by parts as before and applying 1 15. 9ft yields: 

/ (m : (Vv <g> n) ip + B V0 • + v cof W)n - 5 div(§ + z/ cof • n <pj = 0. 

The first term above drops out by (15.111) . Writing Vp = (d T ip)r + (<9^)n we obtain the 
following new boundary equations: 



Nf + u cof *j :{n®n) = 0. 



and: 



d T ((§ + f cof : (n <g> f)j + div ^ + z/ cof * J ra = 0, 
which are, respectively, equivalent to: 

(5.13) § : (n®n) + i/ § : (r <g> r) = onffi, 

(5.14) (1 - i/)9 T : (n ® r)J + div + zv cof n = on dQ. 

In the particular case when (sym n g )2x2 = on dfl, (15.131) and (I5.14p become: 

d 2 m v + u(d 2 TT v-Kd R v)^0 

(2 - u)d T d H d T v + d 3 mH v + K(Av + 2d 2 m v) = 0, 

where K stands for the (scalar) curvature of dQ, so that d T r = Kn. If additionally dQ is 
a polygonal, then the above equations simplify to equations (5) in [TS]. ■ 
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6. Discussion and a proof of Corollary 11.51 

Recall that in classical elasticity it is usually the magnitude of the applied forces or the 
types of boundary conditions which determine the behavior of thin plates or shells, see e.g. 
[TJ. Such exterior constraints are replaced in our case by the geometric constraints induced 
by the prescribed metric. In this line, we conjecture existence of a hierarchy of limit model 
theories, depending on the choice of the tensor a h whose qualifications can be predicted 
by geometric observations. Indeed, in what follows, we demonstrate how a h in (jl.8p is 
related to conditions (11 . 13j) and ( II . 14)) through the Gauss-Codazzi equations or through 
an expansion of the Riemann curvature tensor of the metric G h = (a h ) T a h . In particular, 
these two conditions can be interpreted as the leading order defect in G h from being a flat 
metric. Other choices of scalings in a h should in turn impose the corresponding scalings of 
the energy and the acceptable displacements for the limit model. 

Corollary 11.51 follows now from the next result: 

Lemma 6.1. For any w G W 1 ' 2 ^, R 2 ) and v G W 2 ' 2 (Q,R) , the following are equivalent: 

(i) l g (w,v) = 0, 

(ii) curl((sym /%) 2x2 ) = and curFcurl (e g ) 2x2 = -det((sym /%) 2x2 )- 

The two equations in (ii) are the linearised Gauss- Codazzi-Meinardi equations corresponding 
to the metric Id + 2/i 2 sym(e 9 ) 2x2 and the shape operator /i(sym K ff ) 2x2 on the mid-plate Q. 

Proof. Recall that for a matrix field B G L 2 (Q, R 2 y^) the following two assertions hold true: 

(6.1) curl B = B = -V 2 t> for some v G W 2 ' 2 (Q, R) 



(6.2) curFcurl B = B = sym Vw for some w G W 1 ' 2 ^, R 2 ). 
By ( I6.ip . the first identity in (ii) is equivalent to: 

(sym /%) 2X 2 = -V 2 f . 
Consequently, the second identity in (ii) becomes: 

curFcurl ^ - -Vf <g> Vf + (sym e g ) 2X 2 ) = 0, 
which, in view of (16. 2j) . is equivalent to the existence of w with: 

symVu> + ® ^ v ~ { s J m e g)2x2 = 0. 

Since Q 2 is positive definite on symmetric matrices we see that indeed (ii) is equivalent to 
the vanishing of both terms in X g . 

To identify the equations in (ii), recall the Gauss-Codazzi-Meinardi system [9]: 

d 2 L - dxM = LT\ 2 + M(T 2 12 - r^) - iVr 2 1; 

(6.3) d 2 M - d,N = LY\ 2 + M(r 2 2 - 1*0 - iVr 2 1; 

LN-M 2 = K(EG - F 2 ) 
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which provides the necessary and sufficient conditions for existence of a surface with the 
first and second fundamental forms: 



[da 



PI 



E F 
F G 



and 77 



L M 
M N 



In (16. 3p T l - k denote the Christoffel symbols, and K stands for the Guassian curvature, which 
can be calculated from / and IF Substituting now / = Id+2/i 2 (e / ) 2X 2? II = h(n') 2x2 , where 
e' and k' respectively denote the symmetric parts of e g and K g , and taking into account the 
relations: 

r£ = \g k \d m + d m - d l9ij ) = (^4 + d^ k - d k e' t] )h 2 



0{h 2 



K 



R 



1212 



EG — F 2 

we directly obtain the first identity in (ii): 

h(d 2 {K') n - ^(k')i 2 ) = 0(h 2 ) and h(d 2 {K') 12 - d 1 {K , ) 2 ^j = 0(h 2 ), 
and h 2 det (k / ) 2x2 = -Ri2i2- Now recall that the Riemann curvatures are given by: 

RijM = gim(d k T^ - djTf k + r™r™ - r" fe r™) = d k r i:j - djV ik + o(h A ). 

Hence, after straightforward calcuations, we obtain: 

h 2 det (k') 2x2 = R 1212 = d 2 T 2 u - + 0(/i 4 ) = -/i 2 (curl T curl( e ')2x 2 ) + 0(h 4 
which yields the second identity in (ii). 



7. Appendix A. Approximating low energy deformations: a proof of 

Theorem 11.61 

The first crucial observation follows from the below rigidity estimate, which reproduces 
that of [16] , and it is a non-Euclidean version of the bound in [6] : 

Lemma 7.1. For every u £ W l,2 (Q h , 1R 3 ) and every xq £ Q there exists R £ 50(3) such 
that: 

\ [ \Vu{x) - Ra h (x )\ 2 dx<C (/ h (w) + (diam tt h ) 2 Var 2 (a h )\tt\) . 

The constant C above depends on ||o ft ||ioo ; ||(a /i ) _1 ||£oo, and on the domain Q h . Its depen- 
dence on Q h is uniform for a family of plates which are bilipschitz equivalent with controlled 
Lipschitz constants. 

Proof. Recall that according to the celebrated result in [6], for every v £ W 1,2 (V, M") defined 
on an open, bounded set V C M n , there exists R £ 50(3) such that: 

(7.1) f \Vv-R\ 2 <C V ( dist 2 (Vw,50(3)). 

Jv Jv 

The constant Oy depends only on the domain V and it is uniform for a family of domains 
which are bilipschitz equivalent with controlled Lipschitz constants. 
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In the present setting call A Q = a h (xo) and apply ( 17.1 ft to the vector field v(y) = 
u^Aq 1 !)) G W l ' 2 (Aott h , M 3 ). After change of variables we obtain: 

3ReSO(3) [ \{Vu)A^-R\ 2 <C AoQh [ dist 2 ((Vu)A \ SO (3)). 
Jn h Jn h 

Since the set A Q h is a bilipschitz image of Q h , the constant C Ao Qh has a uniform bound C 
depending on \Aq\, IAq 1 ] and D, h . Further: 



\Vu - RA \ 2 < Cl^ol 4 ^ 1 / dist 2 (Vw, SO(3)A ) 



n h Jn h 



< ClAol^h- 1 (^j dist 2 (Vu(x), SO(3)a h (x)^ dx + j \a h (x) - a h (x )\ 2 dxj 

<C\A^(l*{u) + h- 1 [ 
V Jn 



\a h (x) -a h (x )f 

n h 

The claim follows now through: 

a h (x) - a h (x )\ 2 dx < 2 f (\a h (x) - a h (x')\ 2 + \a h (x') - a h (x )\ 2 ) dx 
(7.2) jQh Jnh 

<C h 2 \d 3 a h \ 2 + ||V tan (o h | n )||| 00 (diain n h ) 2 dx. 



Proof of Theorem 11.61 

The proof follows the line of Theorem 10 [7] (see also Lemma 8.1 [14J). 

1. Let D x i t h = B(x',h) D Q be 2d curvilinear discs in Q of radius h and centered at 
a given x' G Q. On each 3d plate B x i^ = D x i t h x (— h/2,h/2) use Lemma [7.1 1 to obtain 
R x ',h G 5*0(3) such that: 



(7.3) 

< C 



\ [ \Vu h -R x , yh a\x)\ 2 

i^h- 1 j dist 2 {V u h {a h )-\ 50(3)) da; + h 2 Var 2 {a h )\D x , th \J 



with a universal constant C, depending only on the Lipschitz constant of dQ, but indepen- 
dent of h. Notice that we have also used (17. 2p to exchange a h (x') with a h (x) in the left 
hand side above. 

Consider now the family of mollifiers r\ x i : Q — > R, parametrized by x' G Q: 

9(\z'-x'\/h) 



Vx'(z') 



hf Q 9(\y>-x>\/h) dy" 



where 9 G C£°([0, 1)) is a nonnegative cut-off function, equal to a nonzero constant in a 
neighborhood of 0. Then T) x >(z') = for all z' (jL Dx,h and: 

/ T] x > = h' 1 , \\t1x>\\l°° < Ch~ 3 , \\V x >ri x ,\\ L ™ < Ch~ A . 
Jn 
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2. Define now Q h E W 1 ' 2 ^, M 3x3 ): 

Q h ( x >) = [ 7 lx ,{z')Vu h {z)a h {z)- 1 dz. 
By (17. 3p . we obtain the following pointwise estimates, for every x' G f2: 
|QV) - ^,&| 2 < V*>V) {Vu h (z)a h (z)- 1 - R x , ih ) dz) 
( 74 ) <[ \v x >(z')\ 2 dz- [ \Vu h (a h )- l -R x ,, h \ 2 

J V h JB x , h 

[ dM 2 {Vu h {a h )- 1 , 50(3)) dz + h 3 Var 2 (a h )\D x ,, h \ ) 

\VQ h (x')\ 2 = (J^ {V^z')) {Vu h {z)a , \z)- 1 - R xl>h ) dz) 

< [ \V x >vAz')\ 2 dz ■ I \Vu\a h )- 1 - R x , >h \ 2 

J B x',h -' B xi,h 

< Ch- 5 ^ J dist 2 (Vw /l (a fe )- 1 , 50(3)) dz + h z Var 2 {a h )\D xljh \ j . 

Applying the estimates above and in (17. 3p on doubled balls B x >^h we arrive at: 
1 



h 



\Vu h {z)a h {zY 1 - Q h (z')\ 2 dz 



<0(~ f \Vu h {a h )-l-R x ,, 2h \ 2 + ±- f \Q\z') - R x/>2h \ 2 dz J 
(hr x j dist 2 (W> h )-\ 50(3)) dz + h 2 Var 2 (a h )\D x , !2h \j , 



< C 



|VQf < Ch- 2 h~ l / dist'(VM ft (a /l )- i ,50(3)) dz + h 2 Var 2 (a h )\D x ,^ 2h \ 

D x',h \ J B x',2h 

Consider a finite covering £1 = (J -D x / ^ whose intersection number is independent of h (as 
it depends only on the Lipschitz constant of dQ). Sum the above bounds: 

(7.5) i f \Vu h (x) - Q ft (x')a h (ar)| a dx < C (l^(u h ) + h 2 Var 2 (a h )) , 

(7.6) / \VQ h \ 2 < Ch' 2 (I^(u h ) + h 2 Var 2 (a h )) 

Jn 

3. Notice that by (17211 : 
dist 2 (Q/V), 50(3)) < IQV) - < O (h- 2 I^{u h ) + h 2 Var 2 {a h )) ^ as h -)• 0, 
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in view of assumption (4). We may therefore, for small h, project Q h onto SO (3): 

)(3)( 

We further have, by ( 17.5ft : 



i?V)=PsO(3)(QV))- 



/" \R h {x') - Q h {x')\ 2 dx' = [ dist 2 (Q h (x'),SO(3)) dx' 
Jn Jq 



<C7i -1 / \Vu h {x)a h {x)- 1 -Q h {x')\ 2 dx+ / dist 2 {V u h {x)a h {x)' 1 , S0{3)) dx 
\ Jn h Jn h 

< C {Iv(u h ) + h 2 Var 2 (a h )) . 
On the other hand |V-R h | < C\VQ h \ and the claim follows by (1731) and fTTTB]) . ■ 

8. Appendix B: The T-convergence formalism 

Theorems II .31 and II .41 can be summarized using the language of T-convergence [I]. Recall 
that a sequence of functionals T h : X — > R defined on a metric space X, is said to 
T-converge, as h — > 0, to T : X — > R provided that the following two conditions hold: 

(i) For any converging sequence {x h } in X: 

F(limx h ) < liminf F h (x h ). 

(ii) For every x E X, there exists a sequence {x h } converging to x and such that: 

F(x) = \imF h (x h ). 

Corollary 8.1. Define the sequence of functionals: 
T h : W 1 ' 2 ^ 1 ,^) x W^'^R 3 ) x W^'^R 2 ) — > R 

T h (y V,w) = [ ^wivW^a)) */ = fy( x ',t) ~ x ' dt and w = h^Vun, 

+oo otherwise. 



Then T h Y -converge, as h — > 0, to the following functional: 



F(y,V,w) 



X g {w, v) if y(x', t) = x' and V = (0, 0, v) T e W 2 ' 1 
+oo otherwise. 



Consequently, the (global) approximate minimizers of J rh converge to a global minimizer of 
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